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Abstract 

We review some applications of the perturbative technique known as the stochastic limit 
approach to the analysis of the following many-body problems: the fractional quantum 
Hall effect, the relations between the Hepp-Lieb and the Alli-Sewell models (as possible 
models of interaction between matter and radiation), and the open BCS model of low 
temperature superconductivity. 
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I Introduction 



In this paper we review with a certain care the main results concerning three appUcations of 
the so-called stochastic limit approach (SLA), see [1] for the main definitions and some physical 
applications and [2] for some rigorous mathematical results, to three problems in quantum 
mechanics for systems with infinite degrees of freedom. In particular we devote Section II to 
the discussion of the fractional quantum Hall effect, [3]. In Section III we discuss many relations 
between different models of matter interacting with radiation: the Hepp-Lieb and the Dicke- 
Haken-Lax hamiltonian models, and the AUi-Sewell dissipative counterpart, [4]. In Section IV, 
we discuss the open BCS model of superconductivity as originally proposed by Martin and 
Buffet, [5], and some quite recent results related to that, [6, 7]. For reader's convenience, we 
also add a brief Appendix where some crucial facts on the SLA are fisted. 



II The fractional quantum Hall effect 

The Hamiltonian for the quantum Hall effect (QHE) is, see for instance reference [8], 

jj(N) ^ jj(N) ^ ^^^(^) ^ (2.1) 

where Hq^^ is the Hamiltonian for the free N electrons, i?^^^ is the Coulomb interaction: 

fff' = iErr^ (2.2) 

^ i^j \U Lj\ 

and is the interaction of the charges with the positive uniform background. A huge 

literature exists concerning the QHE. We suggest here references [9] and [10]. 
In this section we consider a model defined by an Hamiltonian 

H = Hq^^ + i/o,i? + Aifeb = Ho + XHeb (2.3) 

which is obtained from the Hamiltonian (2.1) by introducing the following approximations: 

(1) the Coulomb background-background interaction is replaced by the free bosons Hamil- 
tonian Hq r given by 

Ho,R^ J io{k)b+{k)b{k)dk (2.4) 

where uj{k) is the dispersion for the free background and b'^{k) and b{k) are bosonic operators 
satisfying the canonical comutation relations: 

[b{k), b+ik!)] = 5{k - U) Hk): b{U)] = [b-^ik): b+{k!)] = (2.5) 
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(2) the Coulomb electron-electron and electron-background interaction is replaced by the 
Prohlich Hamiltonian Hgh, [11], which is only quadratic rather than quartic in the fermionic 
operators: 

Heb^ J ^\r)^{r)F{r- r')(p{'r!)drdr' (2.6) 

where ■ip{r) and 0(r') are respectively the electron and the bosonic fields, while F is a form 
factor. Expanding 0(r) in plane waves, ■0(r) in terms of the eigenstates ipaiz) oi the single 
electron hamiltonian, see (2.14) below, introducing the form factors 

^ (2.7) 

7(2^ 

where 

Vap'ik) := J Jj^)e'^-'-MT)dr (2.8) 
and taking F(r) = e^5(r), [11], we can write 

//e6 = ^ alap{h{gap) + 6+(^)) (2.9) 

a/3 

which is quadratic in the fermionic operators Oq,, a^, 

{fla, 0/3} = {a+, a;^} = {fla, } = 5a/3 (2.10) 
Here we have introduced the smeared operators 

%/?«) = jdkh{k)gpm. (2.11) 

Notice that we are adopting here and in the following the simplifying notation for the 
quantum numbers a := {na,Pa), see formula (2.14) below. 

These are certainly strong approximations. However since from the Frohlich Hamiltonian 
it is possible to recover a quartic interaction, [11], one can say that the Frohlich Hamiltonian 
describes an effective electron-electron interaction which may mimic at least some aspects of 
the original Coulomb interaction. From this point of view it seems natural to conjecture that 
some dynamical phenomena deduced from this Hamiltonian might have implications in the 
study of the QHE. This conjecture is supported by our main result, given by formulae (2.69) 
and (2.70) where we deduce, directly from the dynamics, an obstruction to the presence of a 
non zero x-component of the current, which is quantized according to the values of a finite set 
of rational numbers. This is what we will call a fine tuning condition (FTC). 
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Useless to say that the FTC strongly reminds the rational values of the fiUing factor for 
which the plateaux are observed in the real QHE. We will comment again this fact later. 

In these notes we discuss a model oi N < oo charged interacting particles concentrated 
around a two dimensional layer contained in the {x, y)-plane and subjected to a uniform electric 
field E_ = Ej, along y, and to an uniform magnetic field B = Bk along z. 

The Hamiltonian for the free N electrons Hq^\ is the sum of N contributions: 

TV 

hI,''^^^Ho{i) (2.12) 

i=l 

where Ho{i) describes the minimal coupling of the i-th electrons with the field: 

Ho{^) = ^ {l+-^A{ri)y + eE-u (2.13) 

In the Landau gauge A — —B{y, 0,0), and requiring periodic boundary condition on x, 
'4>{—Lx/2,y) = ijj{Lx/2,y), for almost all y, we find 

V)„p(r) = ^ ^„(y - yf) e^, = huj{n + 1/2) - ^ (eE - (2.14) 

where is the n-th eigenstate of the one-dimensional harmonic oscillator, uj — Ho — 
:^{hku; — eE) and k — ^p, where p e Z. 

Equation (2.14) shows that the wave function ipupiz) factorizes in a x-dependent part, which 
is labelled by the quantum number p, and a part, only depending on y, which is labelled by 
both n and p due to the presence of y^'^ in the argument of the function (/?„. 

Of course the Hamiltonian ifg^^ in (2.12) becames, in a second quantized form, 

^r^ = E£-«X, (2.15) 

a 

where the are the single electron energies, labeled by the pairs a = {n,p) as given in formula 
(2.14) and the are the fermionic operators satisfying (2.10). 

II. 1 The stochastic hmit of the model 

In this subsection we briefiy outline how to apply the stochastic limit procedure to the model 
introduced above. The stochastic limit describes the dominating contribution to the dynamics 
in time scales of the order t/X^, where A is the coupling constant. 
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The starting point is the Hamiltonian (2.3) together with the commutation relations (2.10), 
(2.5). Of course, the Fermi and the Bose operators commute among them. The interaction 
Hamiltonian Heb for this model is given by (2.9) and the free Hamiltonian Hq is given by (2.3), 
(2.4) and (2.15). 

The time evolution of H^h, in the interaction picture is then 

H,b{t) = e^^°*i/e;,e-^^°* = e'J2 <ap{K9o.pe-''^''-'-'^) + b+ig^a^''^^-''"^)) (2-16) 

a/3 

where 

£a(3 = Ea — £l3 (2-17) 

Therefore the Schrodinger equation in the interaction representation is: 

dtUi^^ ^ -iXH,bmi^\ (2-18) 
which becames, after the time rescaling t t/X^, 

dtU^/l = -^^e.(t/A^)^;i. (2.19) 

Its integral form is 

Ui^l. =1- ^ H,b{t'/X')U^//,,dt' (2.20) 
We see that the rescaled Hamiltonian 

J Heb{t/X^) ^e^Y. «'a«^T ^ " ^c.p)gap) + h.c. (2.21) 

a/3 



depends on the rescaled fields 



6a/5,A(t) = ^&(e"^^^"-^«^^^a/3) (2-22) 



The first statement of the stochastic golden rule, [1] , is that the rescaled fields converge (in the 
sense of correlators) to a quantum white noise 

h^p{t) = lim \ b{g^^e-'^ ^^-'-^^) (2.23) 

A— >0 A 

characterized by the following commutation relations 

M), ba'^'it')] = [b+p{t), b+f,,{t')] = (2.24) 

M), bi,^,{t')] = 5e^,,s^,,At - t')G-^'''^' (2.25) 
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where the constants are given by 

Qa^a'p' ^ J^^ J dkg^fs{k)9MSe'<^^-''-''''^ = 2vr / dkg^p{k)gamd{uj{k) - e^p) (2.26) 
The vacuum of the master fields bap{t) will be denoted by rjo: 

bap{t)r}o^O yaP,yt (2.27) 
The limit Hamiltonian is, then, see [3], 

Hit\t) = e'j:{<aMt)+ h-c.) (2.28) 

a/3 

In this sense we say that H^l''\t) is the "stochastic limit" of Hgi,(t) in (2.16). Moreover, the 
wave operator in the stochastic limit satisfies the equation 

dtUt = -iHif\t)Ut (2.29) 

or, in integral form, 

Ut^l-t f H^\t')Ut,dt\ (2.30) 

JQ 

Finally, the stochastic limit of the (Heisenberg) time evolution of any observable X of the 
system is: 

ji(X) = C/+(X(8)lR)C/i, (2.31) 

where 1^ is the identity of the reservoir. Since the hapit) are quantum white noises, equation 
(2.29), and the corresponding differential equation for jt(^), are singular equations and to give 
them a meaning we bring them in normal form. This normally ordered evolution equation is 
called the quantum Langevin equation. Its explicit form is: 

ddtix) = e' Y.{jti[<ap, x]rl^ - r^^[4«- ^])}+ 

a/3 

X]) + MWUp, X])bap{t)} (2.32) 

a/3 



where 



r-^:=E^.«..e„,,,aJa.'G^!^"''^' (2.33) 
Caja'p' ^ £^ J dkg^pik)^:;^^^''^^^-'^-'"'^ = (2.34) 
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_ 1 naPa'P' _ p p / Jk) (k) 1 

2 J - 

The master equation is obtained by taking the mean value of (2.32) in the state 770^"* = 77o®C) 
^ being a generic vector of the system. This gives 



{dtJt{X))^,o = e'j:{M[aUp,X]Tl' - T^'^+fa+a^, X]))^«) (2.35) 



a/3 

and from this we find for the generator 



L{X) = e Y: - <M4^a: (2.36) 

This expressions for L{X) will be the starting point for our successive analysis. 

II. 2 The current operator in second quantization 

The current is proportional to the sum of the velocities of the electrons: 

N J 

JA(t)^a,J2-Ri(t). (2.37) 

i=i 

Here A is the two-dimensional region corresponding to the physical layer, etc is a proportionality 
constant which takes into account the electron charge, the area of the surface of the physical 
device and other physical quantities, and Riit) is the position operator for the i-th electron. 
Moreover N is the number of electrons contained in A. Defining 

N 

XA{t) = J2Riit) , (2.38) 

i=l 

we simply write 

Mt) = ajj,{t) . (2.39) 
Since X\{t) is a sum of single-electron operators its expression in second quantization is given 

by 

XA = ^X^^a+a^ (2.40) 

where 

X^^ = (V^^, Xa^/j^) = I ^l^^{r)ril}^{r)dr (2.41) 

Recall that the ip-yir) arc the single electron wave functions given in (2.14) and that and 
satisfy the anticommutation relations (2.10). In the LLL, = n^ = 0, we find that, [3], 
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= ^S'^^'^P.P. (2-43) 

where y,^,^ := ^(^/S^"^ - y'i'^^) = ^ ^(Pm "Pt)- 

To show how these results can be useful in the computation of the electron current we start 
noticing that, if ^ is a state of the electron system, then 

(/aW), = ac(^^A(i)), = a,{L{Xt,{t)))a = o,cTr{QL{Xt,{t))) (2.44) 

The vector {J\{t))g will be computed now for a particular class of states g, and we will use 
this result to get the expressions for the conductivity and for the resistivity tensor. 

To do this we begin computing the electric current. We first need to find L{X\), L being 
the generator given in (2.36). Since X\ = we have 

L(Xa) =Li(Xa) + /i.c., 

where, as we find after a few computations, 

-^i(-^a) = ^^0,13,^0,' d'^-^"^^ {Xf}ja^a^a'^,aa' - X^aa'^apa'^'aa') (2-45) 

We consider here a situation of zero temperature and we compute the mean value of Li{X^) 
on a Fock A'"-particle state ■0/^ 

iPi = at^... al^^ipo, ik ii,yk I (2.46) 

where / is a set of possible quantum numbers, Nj is the number of electrons in 1 and ■^o is the 
vacuum vector of the fermionic operators, a^V'o — for all a. The order of the elements of / is 
important to fix uniquely the phase of ■^z. Equation (2.44) gives now 

(V^/, JA(t)^/) |i=o= ac(V'/, ^^(^a)^/) (2.47) 

Introducing now the characteristic function of the set /, 

, , f 1 if a e 7 , , 

we get 
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see [3] for the details. Using this equahty, together with 

'^£a/3,£a'a' ~ ^^oc,e/3 <^£„5,£„^/ — <^e/3,£^/ (2.50) 

we find that the average current is proportional to 

{L(Xa))^^ = £i(Xa) + £2(^a) (2.51) 
where we isolate two contributions of different structure: 



>Ci(^a) = E Se.,e,{Xl{0^) - X/(/3)}X/(«')(^a/jG^^"'"' + ^^a^^^"'"'), (2.52) 



£2(Xa) = E Xi{a){l - xiiP')) - G^"' °X/(/3')(l - X/(«))]- 

a/3/3' 



^'/3[G^'''"X7(/5')(1 - X/(«)) - G^--'' Xi{a){i - xm)])- (2-53) 



Using equations (2.42), (2.43) for X» we are able to obtain Ci{Xl') and C2{Xl') for 
i = 1,2. First of all we can show that, even if Ci{X^l^) is not zero, nevertheless it does not 
depend on the electric field. Therefore 

4^i(^f) = (2.54) 

Secondly, the computation of C2{X^l^) gives rise to an interesting phenomenon: due to the 
definition of X^^^ the sum in (2.53) is different from zero only if ^ ppi. Moreover, we also 
must have ep — epr, that is 

27reE 

This equality can be satisfied in two different ways: let us denote TZ the set of all possible 
quotients of the form (n^ — npi)/{ppi —pp)- This set, in principle, coincides with the set of the 
rational numbers. Therefore e 7?.. Then 



np - npi = _ {pD> - pfj) (2.55) 



1) if is not in 7^, (2.55) can be satisfied only ii (3 = (5' . But this condition implies in 

particular that pp = pp', and we know already that whenever this condition holds, then 
X^^^, = 0, so that C2{X'^^) = 0. 
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2) If ^^1^ is in TZ, then we have two possibihties: the first one is again 

which, as we have just shown, does not contribute to £2(-^a The second is 

n^-n^ = (2.56) 

which gives a non trivial contribution to the current. 
Therefore, we can state the following 

Proposition. In the context of our model there exists a set of rational numbers TZ with the 
following property: if the electric and the magnetic fields are such that the quotient 

2TieE 



mw^L 



X 



does not belong to TZ then 

{Ji^\t))^, = 0. 



On the other hand, if condition (2.56) is satisfied, we can conclude that the sum J2a/3i3' ^s/3,s^X- ■ •) 
in (2.53) can be replaced by 

E W(---) = EE'(---) (2.57) 

where means that the sum is extended to all the a and to those (5 and (3' with pp ^ ppi 

satisfying (2.56) (which automatically implies that — Spi). 

Since, as it is easily seen, gai3{k)ga'i3'{M) does not depend on we find that 

where 

^a/3a'/3'^ |° / ci^^.^(^)^;^e^^(^®-«/') (2.59) 

J —oo J 

SO that, using also (2.57), we get 
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where 



-X/(/3')(l - X7(«))(A^"^'" + A^"^'")} (2.61) 

and 

Therefore we conclude that 



im^'X'gl (2.62) 



Let us now compute the second component of the average current: {'^i,L{X\ )tpQ) — 
A(xf) + £2(^f). 

The first contribution is easily shown, from (2.52) and (2.43), to be identically zero, since 

(^6„,6/p„p;3 = ^ap (2.64) 

On the contrary the second term, £2 (-^a ) , is different from zero and it has an interesting 
expression: in fact, due to the factor Sp^^p^, the only non trivial contributions in the sum 
J2pp' ^ei3,spn in (2-53), are exactly those with P = (3'. Taking all this into account, we find that 

>C2(xf ) = Y.i.y^'^ - yS'"^)x/(«)(l - XiimOt^''^ + G^) (2.65) 

aP 

which is different from zero. Furthermore, using (2.58), we get 



were we have defined 



0. = -ppfxi{a){l - xim Im (A^^"^) (2.66) 

a,P 



and A°^"^ is given by (2.59). If we call now 



3x,E- |t=0-«c 
Jy,E- |t=0-«c 
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we obtain the conductivity tensor (see [9]) 



Cf'xx — <^yy — 3y,Ei (^xy — ~(^yx — jx,E (2.67) 

and the resistivity tensor 

'^yy ^xy /q j^qn 

Pxx — Pyy — ~ ~-, Pxy — ~Pyx — o y^-'OO) 

yy ^ ^xy ^yy + ^xy 

After minor computations we conclude that 

r if i n 

Pxy — S mwL^ 6^ -r 27reg ^ 7-) {Z.'0\)) 

^ h ' 2ac<i^Qy mui-'Lx , _ ^ 

1 Qj/ if 27reE ^ -7-) l^- 'Uj 

Let us now comment these results which are consequences of the basic relation (2.56). As it 
is evident from the formula above, the fact that the fine tuning condition (FTC) ( ^^^f G Tl) 
is satisfied implies that p^y ^ 0, so that the resistivity tensor is non-diagonal. Vice-versa, 
if the FTC is not satisfied, then p — Pxx^, 1 being the 2x2 identity matrix. This implies 
that, whenever the FTC holds, then the x-component of the mean value of the density current 
operator is in general different from zero, while it is necessarily zero if the FTC is not satisfied. 

If the physical system is prepared in such a way that e TZ, then an experimental 

device should be able to measure a non zero current along the a;-axis. Otherwise, this current 
should be zero whenever ^TZ. A crucial point is now the determination of the set 71, of 

rational numbers. From a mathematical point of view, all the natural integers Ua and all the 
relative integer p„ arc allowed. However physics restricts the experimentally relevant values to 
a finite set. In fact eigenstates corresponding to high values of Ua and Pa are energetically not 
favoured because the associated eigenenergies e„„p„ in (2.14) increases and the probabilities 
of finding an electron in the corresponding eigenstate decrease (this is a generalization of the 
standard argument which restrict the analysis of the fractional QHE to the first few Landau 
levels). Moreover, high positive values of —pa are not compatible with the fact that Hq must 
be bounded from below, to be a honest Hamiltonian. 

Therefore, in formula (2.56) not all the rational numbers are physically allowed but only 
those compatible with the above constraints. For this reason it is quite reasonable to expect 
that the set TZ consists only of a finite set of rational values. Of course, the determination of 
this set strongly depends on the physics of the experimental setting. 



12 



Finally, let us remark that the sharp values of the magnetic field involved in the FTC may be 
a consequence of the approximation intrinsic in the stochastic limit procedure, which consists 
in taking A — > and t — > oo. In intermediate regions (A > and t < oo), it is not hard to 
imagine that the 5-function giving rise to the FTC becomes a smoother function, and that real 
plateaux, closer to what is observed in the QHE, appear. 

Also, under special assumptions on the 5-dependence of Qx and 9y, together with some 
reasonable physical constraint on the value of the magnetic field, it is not difficult to check that 
Pxx has plateaux corresponding to the zeros of pxy and that, outside of these plateaux, it grows 
linearly with B. 

Ill Laser Models 

In two recent papers, [12, 13], a dissipative laser model has been introduced and analyzed in 
some details. In particular in [12] (AS in the following) the rigorous definition of the unbounded 
generator of the model, which consists of a sum of a free radiation and a free matter generator 
plus a matter-radiation term, is given and the existence of the thermodynamical limit of the 
dynamics of some macroscopic obscrvables is deduced. Moreover, the analysis of this dynamics 
shows that two phase transitions occur in the model, depending on the value of a certain 
pumping strength. In [13] the analysis has been continued paying particular attention to the 
existence of the dynamics of the microscopic observables, which are only those of the matter 
since, in the thermodynamical limit, we proved that the field of the radiation becomes classical. 
Also, the existence of a transient has been proved and an entropy principle has been deduced. 

On the other hand, in a series of papers [14, 15] culminating with the fundamental work 
by Hepp and Lieb [16] (HL in the following) many conservative models of matter interacting 
with radiation were proposed. In particular, in [16] the authors have introduced a model of an 
open system of matter and of a single mode of radiation interacting among themselves and with 
their own (bosonic) reservoirs, but, to simphfy the treatment, they have considered a simplified 
version in which the matter bosonic reservoir is replaced by a fermionic one. This is to avoid 
deahng with unbounded operators. This is what they call the Dicke-Haken-Lax model (DHL 
model in the following). 

In [12, 13] the relation between the AS model and a many mode version of the HL model 
is claimed: of course, since no reservoir appear in the semigroup formulation as given by [12], 
this claim is reasonable but it is not clear the explicit way in which HL should be related to 
AS. Here we will show that the relation between the two models is provided by (a shghtly 
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modified version of) the stochastic hmit (SL). In particular, if we start with the physical AS 
system (radiation and matter) and we introduce in a natural way two reservoirs (one is not 
enough!) for the matter and another reservoir for the radiation, then the SL of the hamiltonian 
for the new system constructed in this way returns back the original AS generator, under 
very reasonable hypotheses. Moreover, the model which we have constructed ad hoc to get 
this generator surprisingly coincides with the HL laser model, [16]. This is the content of 
Subsection 111.1, while in Subsection 111.2 we will consider the SL of the DHL model, [16, 17]. 
We will find that, even if the form of the generator apparently differs from the one by AS, under 
certain conditions on the coefficients which define the model, the equations of motion for the 
observables of the matter-radiation system coincide with the ones given in AS. 

Let us discuss the main characteristics of the three physical models which will be considered 
in this section. In particular, we will only give the definition of the hamiltonians for the HL and 
the DHL models and the expression of the generator for the AS model, without even mentioning 
mathematical details like, for instance, those related to the domain problem intrinsic with all 
these models due to the presence of bosonic operators. We refer to the original papers for these 
and further details which are not relevant in this work. 

We begin with the AS model. 

This model is a dissipative quantum system, S'-^^ consisting of a chain of 2N + 1 identical 
two-level atoms interacting with an n— mode radiation field, n fixed and finite. We build the 
model from its constituent parts starting with the single atom. 

This is assumed to be a two-state atom or spin, Hat- Its algebra of observables, Aat, is that 
of the two-by-two matrices, spanned by the Pauli matrices (ctj,, ctj,, ct^) and the identity, 1. They 
satisfy the relations 

= = = 1; (^x(7y = i(Tz, etc. (3.1) 
We define the spin raising and lowering operators 

= ]j^{crx^icry). (3.2) 

We assume that the atom is coupled to a pump and a sink, and that its dynamics is given by a 
one-parameter semigroup {Tat(t)|tGR+} of completely positive, identity preserving contractions 
of Aat-i whose generator, Lat, is of the following form. 

Lat'7± = -(7iTie)o-±; Lat(Tz = -72((7^ - r)!), (3.3) 

where e(> 0) is the energy difference between the ground and excited states of the atom, and 
the 7's and rj are constants whose values are determined by the atomic coupling to the energy 
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source and sink, and are subject to the restrictions that 



0<72<27i; -l<ri<l. (3.4) 

The matter consists of 2N+1 non- interacting copies of 'Eat, located at the sites r — —N, . .,N 
of the one-dimensional lattice Z. Thus, at each site, r, there is a copy, E^, of E^t, whose algebra 
of observables, Ar, and dynamical semigroup, T^, are isomorphic with Aat and Tat, respectively. 
We denote by (Jr^u the copy of cr^ at r, for u = x,y, z, ±. 

Wc define the algebra of observables, A''^\ and the dynamical semigroup, T^], of the matter 
to be i^^^_j^Ar and ®^_jv^r, respectively. Thus, A^^^ is the algebra of linear transformations 
of C^^"*"^. We identify elements Ar of Ar with those of A^^"^ given by their tensor products 
with the identity operators attached to the remaining sites. Under this identification, the 
commutant, A'r, of Ar is the tensor product ®s^r^s- The same identification will be implicitly 
assumed for the other models. 

It follows from these specifications that the generator, I/^^, of T^] is given by the formula 

ii^i = (3.5) 

where In = {-N, -1, 0, 1, .., N}. Here 

and Lr{ArA'^) = (L^A)A; V^eA, K^A!r (3.6) 

We assume, furthermore, that the radiation field consists of n(< oo) modes, represented 
by creation and destruction operators {a*,ai|/ = 0, . ., n — 1} in a Fock-Hilbert space H-rad as 
defined by the standard specifications that (a) these operators satisfy the CCR, 

[au al] = 5,^1; [a,, aj = 0, (3.7) 

and (b) Tirad contains a (vacuum) vector $, that is annihilated by each of the a's and is cyclic 
w.r.t. the algebra of polynomials in the a*'s. 

The formal generator of the semigroup Trad of the radiation is 

i-rad = Z]"=o^(^^'K^'' ■] + '^i^iat(-)ai - Ki{aiai, .}), (3.8) 

where {.,.} denotes anticommutator, and the frequencies, Ui, and the damping constants, Ki, 
are positive. We refer to [12] for a rigorous definition of Lrad- 
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The composite (finite) system is simply the coupled system, comprising the matter 

and the radiation. We assume that its algebra of observables, B'^^\ is the tensor product 
A'^^^^Tl, where TZ is the *— algebra of polynomials in the a's, a*'s and the Weyl operators. 
Thus, B^^\ like 7^, is an algebra of both bounded and unbounded operators in the Hilbert 
space n^^^ := C^^+^^Hrad- We shall identify elements A, R, of A'^^\ 7^, with A^TLrad and 
^mat'S>R, respectively, with obvious notation. 

We assume that the matter-radiation coupling is dipolar and is given by the interaction 
Hamiltonian 

H^^,^^j:^^Jar,^<f>f^ + h.c.), (3.9) 
where we have introduced the so-called radiation field, whose value at the site r is 

0f ) = -i{2N + l)-V2^;;:^^,a,exp(27riZr/n). (3.10) 

Here the A's are real-valued. A'"— independent coupling constants. 

Among the other results contained in [12] , one of the most relevant is that the map 

is really the generator of a A'"-depending semigroup, T'(^), regardless of the unbounded nature 
of both Lrad and H^^\ 

Let us now introduce the n-modes version of the HL model, [16]. The HL hamiltonian for 
the 2N + 1 atoms and for the n modes of the radiation can be written as follows: 

H^H^^^ + H^^\ (3.11) 

where " S" refers to the system (radiation-|-matter) and " R" to the reservoir. The hamiltonian 
of the system is 

n-l n-1 
n n-l 

+ ^==T E E (^/,+ale-^'^^^-^/- + a,_a,e^-^Vn)^ (3.12) 
V2iV + 1 i^j^ 

[4] . Notice that the presence of P means that we are not restricting our model to the rotating 
wave approximation, (RWA). 
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The hamiltonian for the reservoir contains two main contributions, one related to the two 
reservoirs of the matter and one to the reservoir of the radiation. We have: 

H^R) = H^P)+ Y.H'i''\ (3.13) 
leiN 

where 

n— 1 Tt— 1 

= E / dkuJr,{k)T,{k)^T,{k) + V^^(r](^,)% + r,te)a]), (3.14) 

and 

Hi^^ ^Yl / dkUrn,(k)mli(k)ms^i(k) + ^/a(m\i{hi)ai^- + hx.) + ^/a(m^^ (3.15) 

8 = 1'' ' 

We notice that: 

1) we are using the notation: rj{gj) — J dkrj{k)gj{k) and r^jigj) — J dkrj{k)gj{k). Here dk 
is a shortcut notation for d^. 

2) the functions gj and hi^2 are introduced by HL to regularize the bosonic fields rj{k) and 

^(1,2),/ (^)- 

3) in this model two independent reservoirs, mi;(/c) and m2,;(fc), are introduced for (each 
atom of) the matter, while only one, rj{k), is used for (each mode of) the radiation. This result 
will be recovered also in our approach. 

4) the coupling constanty^ is written explicitly for later convenience. 

The role of each term of the hamiltonian above is evident. We rewrite H as a, sum of a free 
and of an interaction part, in the following way: 

H^Ho + ^Hj, (3.16) 

where 



n-1 



n-1 



Hq^ujrY^ a]aj + cri^z + YY dkum,{k)mli{k)ms,i{k) + Y dkur,j{k)rj{k)^rj{k) 

(3.17) 

and 



n-1 



Hi = E(^](fi)% +rj{9jWj) + Y [("^li(^i)^;,- + h.c.) + {mli{h2)ai^+ + h.c.)] + 
j=o leiN 
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n-1 



n-1 



^ EE (^;,+«]e"'"'''/" + ai^.aje^^'^'/^). (3.18) 



The only non trivial commutation relations, which are different from the ones already given in 
(3.1,3.7), are: 

[rj{k),ri{ky] = 6j,iS{k - k'), [ms,i{k),m\,Jk')] = 5s,s'5i,r5{k - k') (3.19) 



Finally, let us introduce the DHL model. The main difference, [17], consists in the use of a 
fermionic reservoir for the matter, and for this reason the Pauli matrices of both AS and HL 
are replaced by fermionic operators as described in details, for instance, in [17], [4]. Again we 
have 

H^Ho + XHi, (3.20) 

where, this time, 

n— 1 n— 1 . 

^0 = E + E - bl,ib-i) + E y dkujr,j{k)rj{k)^rj{k) + 



I di 



+ E E / dke{k){Bl{k)B,,i{k) + ClmCsM) (3-21) 



and 



j=o 

+ E T.[bUBsA9Bs) + Cs,ii9cs)) + (BligBs) + Cl{9cs))bs,i]. (3.22) 
leiN «=± 

Here gB± and gc± are real function, and the satisfy the following CAR {6±,;, ;} = 1 

and they commutes when localized at different lattice sites: [b±^i, 6± g] — if I ^ s. . 

The commutation rules for the radiation operators (system and reservoir) coincide with the 
ones of the HL model. For what concerns the matter operators (system and reservoirs) the first 
remark is that any two operators localized at different lattice sites commutes, as well as any 
operator of the radiation with any observable of the matter. As for operators locahzed at the 
same lattice site, the only non trivial anticommutators are 

{B^,i{k), Bi^k')} = {C^,i{k), Ci^k')} = 6{k - k'), (3.23) 
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while all the others are zero. Finally, to clarify the different roles between the B and the C 
fields it is enough to consider their action on the ground state of the reservoir (pQ-. 



These equations, together with what has been discussed, for instance, in [17], show that B is 
responsible for the dissipation, while C is the pump. 

III.l AUi-Sewell versus Hepp-Lieb 

We begin this subsection with a pedagogical note on the single-mode single-atom version of 
the AS model. This will be useful in order to show that two reservoirs must be introduced to 
deal conveniently with the matter. After that we will consider the full AS model and wc will 
show that the hamiltonian which produces the AS generator after considering its SL is nothing 
but the HL hamiltonian in the RWA. We will finally comment that adding the counter-rotating 
term (the one proportional to f3 in (3.12)) does not affect this result, since its contribution 
disappear rigorously after the SL. 

The starting point is given by the set of equations (3.3)-(3.10) restricted to n = 1 and 
N — 0, which means only one mode of radiation and a single atom. With this choice the phases 
in ^l"^^ disappear so that the interaction hamiltonian (3.9) reduces to 



and the total generator is L = Lmat + -^^rad + ^[-f^inti •]• 

Let us suppose that the atom is coupled not only to the radiation by means of Hmt, but 
also to a bosonic background m{k) with the easiest possible dipolar interaction: 



(3.24) 



(3.25) 



HMm = (T+m{h) + h.c. 



(3.26) 



Of course this background must have a free dynamics and the natural choice is 




(3.27) 



For what concerns the radiation background the situation is completely analogous: 





H9) + h.c. 



(3.28) 
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are respectively the free hamiltonian and the radiation-reservoir interaction. We take the 
complete hamiltonian as simply the sum of all these contributions, with the coupling constant 
A introduced as below: 

H = Hq + XHj = {fia^ + iUjia'^a + J dk ujm{k)w) {k)m{k) + j dkujr{k)r\k)r{k)} + 
+X{{ar\g) + /i.e.) + {a+m{h) + h.c.) + Ai((7_a^ - /i.e.)}. (3.29) 

Taking the SL of this model simply means, first of all, considering the free evolution of the 
interaction hamiltonian, Hi{t) = e*^°*i7/e~*^°*. It is a simple computation to obtain that, if 
cur = 2/1, then 

Hi{t) = (ar^(^e^('^'^ -'"«)*) + h.c.) + ((7+m(/ie^(2'*-'^'")*) + h.c.) + Ai((7_a^ - h.c). (3.30) 

In this case the SL produces, as discussed in detail in [4], the following effective time-depending 
interaction hamiltonian: 

Hi'^\t) = (arj(t) + h.c.) + {a+mhit) + h.c.) + i{a-a^ - h.c), (3.31) 

where the dependence on A disappears and the operators rg{t), mh{t) and their hermitian 
conjugates satisfy the following commutation relations for t > t', 

hit),rlit')] = r^f^sit - 1'), K(t),mi(t')] = r^-^sit - 1'). (3.32) 

Here we have defined the following complex quantities: 

T^3) ^ I dk\g{k)\^e-"-^^'-''^-^^^\ t'l^ = f dr j dk\h{k)\^e-'^^^-^^^^^^\ (3.33) 

We want to stress that the restriction t > t' does not prevent to deduce the commutation rules 
(3.36) below, which are the main ingredient to compute the SL. However, the extension tot <t' 
can be easily obtained as discussed in [1]. Of course the functions h and g must be chosen in a 
such way that the integrals above exist finite! 

In order to obtain the generator of the model we introduce the wave operator Ut (in the 
interaction representation) which satisfies the following operator differential equation: 

dtUt^ -iHf\t)Uu ^iih C/q = 1- (3.34) 

We have already commented in Section II that, at least for a large class of quantum mechanical 
models, the equation above can be obtained as a suitable limit of differential equations for a 
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A-depending wave operator, [1]. Analogously, rg{t) and mh{t) can be considered as the limit (in 
the sense of the correlators) of the rescaled operators ^r{ge"^'''^''~'^"^^/^^) and im(/ie*(^'*~'^"')*/'*'^) . 
It is not surprising, therefore, that not only the operators but also the vectors of the Hilbert 
space of the theory are affected by the limiting procedure A — > 0. In particular, the vacuum 770 
for the operators and m/j, mh{t)r]o — rg(t)r]Q — 0, does not coincide with the vacuum ipo for 
m{k) and r{k), r(k)(po = m{k)(po = 0, see [1] for more details. 

Equation (3.34) above can be rewritten in the more convenient form 

Ut^TL-i t Hf(t')Ut'dt', (3.35) 
Jo 

which is used, together with the time consecutive principle, [1], to obtain the following useful 
commutation rules 

[rg{t), Ut] = -iT^f^aUt, K(t), Ut] = -irL'^^.C/*. (3.36) 

If we define the fiow of a given observabe X of the system as Jt(X) = Ul{X ® iRjUt, the 
generator is simply obtained by considering the expectation value of dtjt{X) on a vector state 
r/o^^ = riQ ^ ^, where ^ is a generic state of the system. Using formulas (3.34), (3.36) and their 
hermitian conjugates, together with the properties of the vacuum r]Q, the expression for the 
generator follows by identifying L in the equation {dtjt{X)) (5) = {jt{L{Xj)) n). The result is 

L{X)^L,{X) + L2{X) + Ls{X), 

Li{X) = r^f\a^,X]a-T^_!^a^[a,X], L^iX) = rL^)[(7+, X](7_ - rL^V+[(7_, X], 

L^{X) = i^la^a) - cr+a, X] (3.37) 

It is evident that both Li and L3 can be rewritten in the same form of the radiation and 
interaction terms of the AS generator but this is not so, in general, for L2 which has the form 
of the AS matter generator only if the pumping parameter 77 is equal to —1. 

This is not very satisfactory and, how we will show in the following, is a consequence of 
having introduced a single reservoir for the atom. We will show that the existence of a second 
reservoir allows for the removal of the constraint r] — —1 above. 

With all of this in mind it is not difficult to produce an hamiltonian which should produce 
the full AS generator for the physical system with 2A'" + 1 atoms and n modes of radiation. 
With respect to the one discussed above, it is enough to double the number of reservoirs for the 
matter and to sum over I e 1^ for the matter and over j — 0, 1, n — 1 for the modes. The 
resulting hamiltonian is therefore necessarely very close to the HL one: 

H = Ho + XHi, (3.38) 
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with 

n— 1 2 „ n— 1 „ 

(3.39) 

and 

n-1 

Hi = +T'ji9j)a]) + [(^L(^i)^i - + ^-c-) + {■m{i(h2)cri,+ + h.c.)] + 

3=0 I&In 

+A^((/.J^V,,+ + /i.c.), (3.40) 

where the radiation field has been introduced in (3.10). It is clear that, but for the RWA 
which we are assuming here, there are not many other differences between this hamiltonian 
and the one in (3.11)-(3.15). It is worth mentioning that A appears both as an overall coupling 
constant, see (3.38), and as a multiplying factor of I]zG/jv(*i^/^'''^^+ ^ h.c.) and plays the same 
role as ^/a in the HL hamiltonian. As for the commutation rules they are quite natural: but 
for the spin operators, which satisfy their own algebra, all the others operators satisfy the CCR 
and commute whenever they refer to different subsystems. In particular, for instance, all the 
m\ i{k) commute with all the m\^i,{k')., for all A;, k' and Z, V . 

The procedure to obtain the generator is the same as before: we first compute Hi{t) — 
^iHotjj^^-iHot^ which enters in the differential equation for the wave operator. Taking the limit 
A — > of the mean value in the vector state defined by ifl^^ = v'o ® of the first non trivial 
approximation of the rescaled version of Uf we deduce the form of an effective hamiltonian, 
H^^''\t), which is simply 

Hi''\t) = E {cTi,^m,,i{t)+h.c.)+ E {ai,^m,M+h-c.)+ E {<pf^cTi,^+h.c.). 

j=o leiN leiN leiN 

(3.41) 

Again, we are assuming that cur = 2/i, which is crucial in order not to have a time dependence 
in the last term of Hi'''\t) in (3.41). 

The only non trivial commutation rules for t > t' for the new operators are: 

[r.jW,4'(0] = rL^k/'^(t-0, 

[m,,iit),mli,it')] = r^!!'^5i,i,5it - t'), (3.42) 

[m,,i{t),my{t')] = T^^'^6i,,6{t-t'), 



22 



where we have defined the following complex quantities: 

J —OO J J —oo J 

rL'^^) = r dr f rf/t|/i2(/t)pe^('"'-2(fe)+'-HK. (3.43) 

J—oo J 

Prom the commutation rules above and since Ut^l-iJ Hf\t')Ut>dt' , we get 
VgAi)^ Ut] = [m,,i{t), Ut] = -zr^Vz,_C/t, [m2,i{t), Ut] = -iV^-''^ai,+Ut. (3.44) 

The expression for the generator can be obtained as for the N — n = \ model described 
before, that is computing the mean value (dtjtiX)) (e). Here, as before, jtiX) is the fiux of the 
system observable X, jt{X) = Ul{X ® Ir)^/*, and 770 is the vacuum of the operators rgj{t) and 
uis^iit), s = 1,2. The computation gives the following result, which slightly generalize the one 
in (3.37): 

L{X)^L^{X) + L2iX) + L-,{X), 

n-l 
j=0 

L,{X) = t E [(^S'^V+.i + h.c.),X]. (3.45) 

It is not difficult to compare this generator with the one proposed by AS, see formulas ((3.3)- 
(3.10)), and the conclusion is that the two generators are exactly the same provided that the 
following equalities are satisfied: 

^r^9).^^.^ sRrL^). = fc,-, K(rL'^^^ + rL'''^) = 71, Q{r^y - r^!!'^) ^ e 

sRr^) = 1^2(1 -r?), ^r^^^^^^^,{i + n). (3.46) 

Therefore, if we start with the HL hamiltonian, choosing the regularizing functions in such a 
way that the equalities (3.46) are satisfied, the SL produces a generator of the model which is 
exactly the one proposed in [12, 13], with the extra minor constraint 72 = 271, which is a direct 
consequence of (3.46). 

We conclude this subsection with a remark concerning the role of the RWA and its relation 
with the SL. In particular, this is a very good approximation after the SL is taken. To show 
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why, we first notice that adding a counter-rotating term (extending the one in (3.18)) to the 
interaction hamiltonian Hj in (3.40), considering the same couphng constant for both the rotat- 
ing and the counter- rotating term {(3 — a), simply means to add to Hi in (3.40) a contribution 
hke XY.ieiNi^i^^^i- + h.c). While the rotating term, if 2/x = our, does not evolve freely, the 
free time evolution of this other term is not trivial. However, the differences with respect to the 
previous situation all disappear rigorously after the SL, because these extra contributions to 
the mean value of the wave operator go to zero when A — > 0, so that at the end the expression 
for the generator is unchanged. This allows us to conclude that the full HL hamiltonian is 
equivalent to the AS generator, where the equivalence relation is provided by the SL. 

III.2 The SL of the DHL model 

We now consider the SL of the DHL model and we will get the expression of the related generator 
showing that, under some conditions on the quantities defining the model, the equations of 
motion do not differ from the ones in AS. The free evolved interaction hamiltonian Hj in (3.22) 
is, 

n— 1 
j=0 

+ E [bUB+,ii9B+e''^'-'^) + C^,i{9c+e'"~'-'^)) + (^^.(^B+e-*^'^-)) + Ci,,(5c+e-*(--)))6+,; + 
leiN 

+bl,iB_,iigB-e-''^^+'^) + C_,iigc-e-''^^+'^)) + (Bl.igB-e''^^^'^) + Cl,igc-e''^^'-'^))b_,i]. (3.47) 

Following the usual strategy wc conclude that (the rescalcd version of) the wave operator 
Ux{t) = 1 — iX H i{t')U\{t')dt' converges for A — >• to another operator, which we still call a 
wave operator, satisfying the equation 

Ut^l-i / Hf'\t')Ut'dt', or, equivalently dtUt = -iHf'\t)Ut, with Uq = 1. (3.48) 

Jo 

is an effective time dependent hamiltonian, found using the usual strategy, defined 

as 

n-l 

Ht\t) =E(%4W + h-c.) + + h.c.) + E [bUf^+,iit) + 7+/^)) + 

j=o leiN 

HPUt) + ^Ut))K,i + bliW-At) + 7-,/ W) + iPUt) + (3-49) 

The operators of the reservoir which appear in are the stochastic limit of the original 

(rescaled) time evoluted operators of the reservoir and satisfy (anti-)commutation relations 
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which are related to those of the original ones. In particular, after the SL, any two operators of 
the matter (system and reservoirs) localized at different lattice sites commutes, as well as any 
operator of the radiation with any observable of the matter. As for operators localized at the 
same lattice site, the only non trivial anticommutators are 

{M).Pliit')} = ^(t - t')^-^"^^^ b±At):7Ut')} = - (3-50) 
which should be added to 

KAt),rgAt')] = kfS(t - t')r%. (3.51) 

In all these formulas the time ordering t > t' has to be understood and the following quantities 
are defined: 



' J —oo J J —oo J 

= r dr [ dk{gc±{k)fe-'^'^''^^^'^\ (3.52) 

J —oo J 

We call now 770 the vacuum of these limiting operators. We have 

r,,,(t)r/o = l3±At)rio = ll,iit)rio = 0- (3.53) 

Paying a little attention to the fact that here commutators and anti-commutators simultane- 
ously appear, we can compute the commutators between the operators rgj(t)'*, 7j_ ^(t), /S^^ i(^) 
and the wave operator Ut by making use of (3.50,3.51). We give here only those commutation 
rules which are used in the computation of the generator: 

KAt)^ Ut] = -iV^l]^a^Ut, [Mt), Ut] = -iV^^''h^,iUt, [ii^i{t), Ut] = iV^-^%iUt. (3.54) 

Using the usual strategy, and restricting to quadratic matter operators for technical reasons 
(this condition can be avoided), [4], we get 

L(X) = Li(X) + L2(X) + L3(X), 

3=0 

L2{x) = Y.{^^^''\bl,uX]h^,-v^^'-hl,[h^^ + 

L,{X) = I Y: [{<l>f'%ih-,i + h.c.),X]. (3.55) 
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We see that the first and the last terms exactly coincide with the analogous contributions of the 
AS generator, but for a purely formal difference which is due to the different matter variables 
which are used in the two models. The second contribution, on the other hand, cannot be 
easily compared with the free AS matter generator. What is convenient, and sufficient, to get 
full insight about L2, is to compute its action on a basis of the local algebra, that is on 
(= (7+,/) and on b\_ fi+^i — bl_ fi-^ (= (Tz,i), all the others being trivial or an easy consequence of 
these ones. It is not hard to find the result: 

_ic^[r(_^+) - rL^-) - rL^+) + rL^-^]), 

L2{bl^fi+,i - blib^,i) = 2{-bl^fi^,i{^T^^^^ + rf +)) + KrL^+) + 

+bl^ib.,i{^[T^^-^ + rL^-^) - 3firL^-^). (3.56) 

The equation for a+^z is recovered without any problem, modulo some identifications (3?[rL^+^ + 
r^"-* + r^'*''' + r^"-*] = 71, ...), while to recover the equation for az,i it is necessary to choose 
properly the regularizing functions which define the different r_. In particular we need to have 
the following equahty fulfilled: 

sjfj(r(s+) ^ p(c+).| _ s)fj(rL^-) + rL^-)). (3.57) 

Under this condition we can conclude that the SL of the DHL model produces the same differ- 
ential equations as the AS generator, as already happened for the HL model. It is also easy to 
check that, as a consequence of our approach, we must have 71 = 72 in the generator we obtain. 
Of course this result is not surprising since already in the HL paper, [16], the fact that the two 
models are quite close (under some aspects) was pointed out. Here we have learned also that 
the SL of both these models, at least under some conditions, give rise to the same dynamical 
behavior. 

IV The open BCS model 

In this section we review some results on the Open BCS model, originally introduced in [5, 17], 
obtained using the SLA, [6]. 

The main outcome is that the same values of the critical temperature and of the order 
parameters can be found using the SLA, in a significantly simpler way. This simplification 
allows us to focus our attention on some aspects of the model which could appear not so clearly 
using the original approach. 
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IV. 1 The Physicals Model and its stochastic limit 

Our model consists of two main ingredients, the system, which is described by spin variables, 
and the reservoir, which is given in terms of bosonic operators. We refer to Section III for 
the definition of the relevant operator algebras. The system is contained in a box of volume 
V — L^, with N lattice sites. We define, following [5, 17] 

N N 

H'r'-e^-^-j.j:4-7, (4.1) 

j=l i,j=l 

where the indexes i, j represent the discrete values of the momentum that an electron in a fixed 
volume can have, o"^ creates a Cooper pair with given momentum while a~ annihilates the 
same pair, e is the energy of a single electron and —g < is the interaction close to the Fermi 
surface. As we can see, only the ± component of the spin, that is the x, y components, have 
a mean field interaction, while the z component interacts with a constant external magnetic 
field. The algebra of the Pauh matrices is given by 

= Si^al [a^, aj] = T2S,,a^. (4.2) 

We will use the following reahzation of these matrices: 

If we now define the following operators, 

1 ^ 

^N-i^i: <, Rn = St,S], = Ri, (4.3) 
-'^ i=i 

if^^*^ can be simply written as if^^*^ = N{eS%—gRN) and it is easy to check that the following 
commutation rules hold: 

[S%,Rn] = [Hn^''\Rn] = [Hn^''\S%] = 0, 

for any given N > 0. It is also worth noticing that the intensive operators are all bounded 
by 1 in the operator norm, and that the commutators [S^, ctj] go to zero in norm as when 
N — > 00, for all j, a and /3. 

Our construction of the reservoir follows the same steps given in [17], but for the commu- 
tation rules. We introduce here as many bosonic modes a^j as lattice sites are present in V. 
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This means that j — 1,2, ...,N. pis the value of the momentum of the j-th boson which, if we 
impose periodic boundary condition on the wave functions, has necessarily the form p — ^n, 
where n — {rii, n2, n^) with rij e Z. These operators satisfy the following CCR, 

[af,i, a^,j] = [at, at] = 0, [a^-^, at] = SijS^^ (4.4) 

and their free dynamics is given by 

H^r^ = E E 4«P.- (4-5) 

Here A^v is the set of values which p may take, according to the previous remark: — {p — 
^n, n e Z!^}. It is useful to stress that the energy of the different bosons is clearly independent 

of the lattice site: ^p — ^ — ^^2ml^^^^^ ■ 
The interaction is given by 

^5.'^ = E«%(/) + ^-c.), (4.6) 

i=i 

where we have introduced aj{f) — J2peAN '^pjfip)^ f being a given test function which will be 
asked to satisfy some extra conditions, see equation (4.24) below and the related discussion. 
The finite volume open system is now described by the following hamiltonian, 

Hj, = H% + XH^^\ where K = ^J^'^^ + H^,''^ (4-7) 

and A is the coupling constant. 

As usual, the first step in the SLA is the computation of the free evolution of the interaction 
hamiltonian: 

Hii\t) = e^^^*i/S;)e-^^* = E(e^''-^^V;e-<''^'*e<^^\,(/)e-<^^^* + h.c). (4.8) 

The computation of the part of the reservoir is trivial and produces 

e'^N '*a,(/)e-'< = a,(/e-**^), 

where aj{fe~^^^) — JZpeApf This is an easy consequence of the CCR (4.4). The 

free evolution of the spin operators is more difficult and its expression can be found in [5, 17], for 
instance, where it is shown how to obtain the time evolution in a semidassical approximation. 
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i.e., when the free time evolution of the intensive operators S% are replaced by their limits in 
a suitable topology, [19] and [6]. 

The differential equations of motion for the spin variables are 

^ = 2tgial{t)S],{t)-ajit)SUt)). 
where a^{t) = g^^'V^e"^^^"^'* and S%{t) = e^^^'^'^-S-^e"^^^"^'* = ^ Ef=i e^^^'^'V^^e"^^^"^'* = 

Let us now call 8°" — !F — strong \imf^^^ S^, where ^ is a suitable family of vectors. 
The proof of the existence of this limit (together with all its powers) may be found in [19] and 
references therein. We can now take the sum over j = 1, 2, N of (both sides of) the equations 
in (4.9), divide the result by N, and then consider the — stronglimi^^^ of the equations 
obtained in this way. We find that S^{t) = and S+{t) =i{2e + gS^{t))S+{t). These equations 
can be easily solved: S%t) = 5° = {Sy and S+{t) = ^+e^(2^'+s5")*. Of course S-{t) = {S+{t)y. 
The system (4.9) becomes, if we replace Sf^{t) with its — strong limit S°'{t), 

'4^ = 2^ia|{t) + ^gS^t)a^{t) 
'4^ = 2^g{a|{t)S-{t)-a-{t)S^t)). 

This system is called the semiclassical approximation of (4.9), and it can be explicitly solved: 

a+(t) = e'^V'o + e'^'+'^^Vi + e'^""^^' p>_, (4.11) 

where we have defined the following operators 

Pi = 4? (2-5-^+ + + 2S^a-) 

P% = ^ {aS-^aj + - gS-a-) (4.12) 

and the following quantities 



CO = g^{Sy + iS+S-, u = 2i+ gS'^. (4.13) 

Defining further 

I'aip) ^ - + aou, (4.14) 
where a takes the values 0, + and — , the operator Hj^\t) in (4.8) becomes 

N 

Hii\t) = E E {p'a^Afe''''') + h.c) . (4.15) 

j=l a=0,± 
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The next step in the SLA consists in computing the following quantity 

hit) = jyt, dhujtot , (4.16) 

and its limit for A going to zero. Here the state ujtot is the following product state cutot — ^sys ^p-i 
where ujsys is a generic state of the system, while ujp is a state of the reservoir, which we will 
take to be a KMS state corresponding to an inverse temperature (3 = It is convenient 
here to use the so-called canonical representation of thermal states, [1], which is sketched in 
the Appendix. Then we introduce two sets of mutually commuting bosonic operators {c^j}, 
7 = a, 6, as follows: 

ap, = V^^S + V^^S'^ (4-17) 

where 

J. 1 i e-/3ef? 

m{p) = ujp{ap,jal>) = ^ _ n{p) = ujp{al^ap,j) = (4.18) 

The operators c^°- satisfy the following commutation rules 

[cg,cg']=5,,5,Va., (4.19) 

while all the other commutators are trivial. Furthermore, we introduce the vacuum of the 
operators c^!*, $o- 

cg$o = 0, Vp e A^, 3 = 1, ..N, a^a,b. (4.20) 



Finally if we define fmip) = \Jm{p)f{p) and fn{p) = \/n{p)f{p),we get 

a,.(/e^*^«) = cf\fme'"'") + cf\u^'^-), (4.21) 
using the usual notation for cf\g) and cf^\g). Therefore we have 

Hii\t) = E E {pi (c^Hfme''''") + Vne^*^i) + h.c] , (4.22) 

and the KMS state LOp can be represented as the following vector state, as in a GNS-like 
representation: 

u;fs(Xr)^{^o,Xr^o), (4.23) 

for any observable of the reservoir, X^, since LOf^ is a gaussian state, [1]. This fact, together 
with (4.20) and with the commutation rules (4.19), simplifies the computation of the two point 
function Utot (^H^j^\^)Hj^\^)^ , which produces 

u^tot(Hi^\hH^^\'^,)) = j: E E H..(p^y;)i/-(piiv^-(^"^e-^''^(-v 

^ ^ ^ ^ i=la,/3=0,±peAjv 
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Since we are interested to the limit A — of Ix{t) we need to impose some conditions on the 
test function f{p), [1]. In particular, we will require that the following integral exists finite: 



/ dr J2 Ifrm'^^'^""^^ < oc, (4.24) 
where fr{p) is fm{p) or fn{p) and I'aip) is given in (4.14). Under this assumption we find that 

N 

I{t) = lhn/,(t) = -tE E {^sys{f^af^a^)ri-^ +UJsys{f^afi)r^^} , (4-25) 



j=l a=0,± 

where the two complex quantities 



rL"^= r E l/m(p)re-'^'^"®, r dr E lUm^e''""^^^ (4.26) 



both exist because of the assumption (4.24). 

To this same result we could also arrive starting with the following stochastic limit hamil- 
tonian 

Hi^'\t) = E E {/^« tew + c^fit)) + h.c] , (4.27) 

j = la=0,± ^ ^ ^ ' 

where the operators (t) satisfy the following commutation rule, 

[c%ht), 4f{t')] = Sj, S,,S{t - t')rW, for t > t'. (4.28) 

We mean that, as it is easily checked, the following quantity 

J{t) - i-tf fdt, r dt2ntot{H^f{h)H^N\t2)) 

Jo Jo 



coincides with I{t). Here ^2^* = ^sys ^ = ^sys (^o, ^o), where is the vacuum of the operators 

cj^t): ci].^(t)^o = for all a, j,7 and t, [1]. 

We now use H^^^\t) to compute the generator of the theory. Let X be an observable of 
the system and 1^ the identity of the reservoir. Its time evolution (after the stochastic limit is 
taken) is jt{X ® 1^) = Ul{X ® lr)^t; where Ut is the wave operator satisfying the following 
differential equation dtUt = —iH^''\t)Ut, whose adjoint is dtU} = iU}H^''\t). 

Then we find 

dtJtiX ® 1,) = tU^[Hi^'\t),X(^ IrPt = 
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j=la=0,± ^ 



Using now the commutation rules 



:]{t),u,] = -I jy:]{t),H^\t')]u,dt' = -I fyjT^:^5{t-t'))u,dt' = -ip>Jt^:^u, (4.29) 



and 

and their adjoints, we find that 

N 



[c^ah)M--ip'a^^a^Ut, (4.30) 



j=l a=0± ^ ^ 



which has to be computed on the state '^tot- Therefore, since the generator L satisfies the 
equahty ntot{dtjt{X (g) 1^)) = Q.tot{jt{L{X))), we get 



L{x) = E E xyJ^^:^ + k\ x]pirw - xjri") - x]ri'') (4.3i) 

This expression can be made simpler if the observable X is self-adjoint [X = X^) . In this case 
we have 

L{X) = L^{X) + L2{X), (4.32) 

where 

N N 

^iW = E E + L,{X)^Y.T.H\xW'a^ + h.c]. (4.33) 

j=l a=0± j=l a=0± 

IV. 2 The phase transition 

As discussed in [5, 17], S% and are the relevant variables whose dynamics must be considered 
to analyze the phase structure of the model. These intensive operators arc both sclf-adjoint, 
so that we can use equations (4.32) and (4.33) instead of (4.31). As a matter of fact, in both 
[5] and [17] these equations of motion are considered only as an intermediate step to compute 
the equation for A^v = \^n^ i which is called the gap operator. We have shown in [6] that the 
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same conclusions as in [5, 17] can be obtained without introducing A^v but working directly 
with Rn and S%. 

First we focus on L(S%) = Li(S%) + L2(S%). We have, using (4.2), (4.3) and (4.33) 

j=l j=l a=0,± 

whose limit in the T — strong topology exists, [19] , and is given by 

T rQQ\ T + V T (Q^\ 8g^'g°('g+'g~)^ [a>T^(a) , ^n(a) ^ + ^ ] 

U{S ) T -strong \.^^U{S^) = + ^^J3^| ' 

(4.34) 

where 3?rf^^ indicates the real part of , [6] . 

The computation of L2{S^) :— T — strong\\m.N-^oo ^2(8^) follows essentially the same steps 
and produces 



[6] , so that the final result is 



Here we have introduced, for brevity, the function 



(a) , ^Ua) ^ + 9 , sorW ^ + 9 , ^Ub) ^-9 



h{s\ S+S-) - kf; ^^^^2 + ^r- _ ^^0)2 + ^^^^^0)2 + ^r. _ , 

(4.37) 

and we have written explicitly the dependence of h on S^S^ = T — strong limN^oo S^S^, 
see (4.13). It is interesting to observe that the same function h{S^,S^S~) appears in the 
computation of Li^S^S") := T — strong limN^oo L{S^S]:f). Again, since (S'^S'^)''" = S^S^, 
we can use formulas (4.32) and (4.33). Here the computations are significantly harder, but no 
difficulty of principle arises. As a technical tool it is convenient to use the fact that, in the 
limit A?" ^ oo, all the intensive operators commute with all the local operators of the system, 
limjv^oo['S'^; Cj] = 0, for all a, (3 and j. Therefore we get 

L{S+S-) = -^-^^^^^^-Lh{S\ S+S-). (4.38) 

The phase structure of the model is now given by the right-hand sides of equations (4.36) 
and (4.38), see [5, 17], and, in particular, from the zeros of the functions 

h{x, y) = -^^-^h{x, y), /2(x, y) = -^^^-^h{x, y), (4.39) 
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where we have introduced, to simphfy the notation, x — S'^ and y — S'^S~. In particular, the 
existence of a superconducting phase corresponds to the existence of a non trivial zero of /i 
and /2, [5, 17]. Due to the definition of /i and /2 it is clear that a pair {xo,yo), with Xo ^ 
and yo ^ 0, is such that fi{xo, yo) — f2{xo, yo) = if and only if it is a zero of the function h: 
h{xo,yo) = 0. In order to find such a zero, it is first necessary to obtain an explicit expression 
for the coefficients 3fJr±^. This is easily done using the definitions in (4.26), since we get 



1 POO 

^ / E \Up)\'e-''^^^ = TT ^ lUm'SiMp)), (4.40) 



and 



^rf^n E \fn{p}miy±{p})- (4.41) 

It is now almost straightforward to recover the results of [5, 17]. Following Buffet and Martin's 
original idea, we look for solutions corresponding to u — 0. This means that, because of (4.13), 
the value oi x — is fixed: x — —2e/g. Moreover, with this choice, iy+{p) — cu — e^, which 
is zero if and only if a; = ep. Also, we have iy-{p) — —cu — ep, which is never zero. For these 
reasons we deduce that 3?rL^^ = 0, J — a,b, while the sums in (4.40) and (4.41) for 3?r+ are 
restricted to the smaller set, C Ajv, of those values of p such that, ii q E £n then — uj. 
Therefore, recaUing the expression of m{p) and n{p) in (4.18), we find 

- E \m\\ ^r? = TT^^ E \m?^ (4.42) 



so that equation /(x, y) looks like 

= 

or 



= 9±^_ (4.43) 
g-uj 



I3lu 



LO. 



This equation is the crucial one, which replaces the one obtained in [5, 17], gftanh 
As a matter of fact, in [6] it is also proven that these two equations are equivalent, and for that 
we recover exactly the same values of the critical temperature and of the order parameters as 
in [5]. 
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IV.3 More results 



We have shown how the SLA can be successfully used to analyze the phase structure of low 
temperature superconductivity analyzing a strong coupling BCS model, considered as an open 
system interacting with a bosonic thermal bath. 

This procedure is rather direct and technically much simpler than the one used in [5]. 
Among the other simplifications, for instance, a single equation h{x, y) — must be solved 
instead of the system fi{x, y) = f2{x, y) = 0. This suggests to use the SLA also to modify the 
original model in the attempt of getting some insight on high-temperature superconductivity. 
This project started quite recently, [7] , by introducing two reservoirs instead of only one, as we 
did here, to see whether the value of the critical temperature increases because of the presence 
of this second reservoir. Our results seem rather promising but not jet definitive. A deeper 
analysis is being presently undertaken. 
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Appendix: Few results on the stochastic Hmit 



In this Appendix we will briefly summarize some of the basic facts and properties concerning 
the SLA which are used all throughout the paper. We refer to [1] and references therein for 
more details. 

Given an open system S + TZ we write its hamiltonian H as the sum of two contributions, 
the free part Hq and the interaction XHj. Here A is a coupling constant, Hq contains the free 
evolution of both the system S and the reservoir TZ, while Hj contains the interaction between 
S and TZ. Working in the interaction picture, we define Hj[t) = e*^"*if/e~*^'^* and the so called 
wave operator Ux{t) which is the solution of the following differential equation 

dtUxit) = -iXHi{t)Uxit), (A.l) 

with the initial condition Ux{0) = 1. Using the van-Hove rescaling t — > see [17, 1] for 
instance, we can rewrite the same equation in a form which is more convenient for our pertur- 
bative approach, that is 

dtUxi^,) = -jHiiy^Wxiy,), (A.2) 
with the same initial condition as before. Its integral counterpart is 

which is the starting point for a perturbative expansion, which works in the following way. 

Suppose, to begin with, that we are interested to the zero temperature situation. Then let 
ipo be the ground vector of the reservoir and ^ a generic vector of the system. Now we put 
•^0^^ = V^o ® We want to compute the limit, for A going to 0, of the first non trivial order of 
the mean value of the perturbative expansion oiU\{t/)^) above in that is the limit of 

h{t) = j^dt^ C MHi{%)HA%))^^,, (AA) 

for A 0. Under some regularity conditions on the functions which are used to smear out the 
(typically) bosonic fields of the reservoir, this limit is shown to exist for many relevant physical 
models, sec [1], and [3, 4, 6] for few recent applications to quantum many body theory. It is at 
this stage that all the complex quantities like the F^'^'^'s we have introduced in the main body 
of this paper appear. We define I{t) = lim;v^o -^a (^) • In the same sense of the convergence of 
the (rescaled) wave operator U\{^) (the convergence in the sense of correlators), it is possible 
to check that also the (rescaled) reservoir operators converge and define new operators which 
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do not satisfy canonical commutation relations but a modified version of these. For instance, 

— * 

in Section II this procedure has produced the operators ha^p{t) starting from h{k). Moreover, 
these limiting operators depend explicitly on time and they live in a Hilbert space which is 
different from the original one. In particular, they annihilate a vacuum vector, 770, which is no 
longer the original one, ip^. This is what happens, for instance, if </7o depends on A, </7o ^o'K 
and considering r/o as the following limit: 770 = limA^o V^o'*'^ • 

It is not difficult to deduce the form of a time dependent self-adjoint operator Hj^''\t), which 
depends on the system operators and on the limiting operators of the reservoir, such that the 
first non trivial order of the mean value of the expansion of f/f = 1 — i Jq Hj^''\t')Ut'dt' on the 
state rjQ^^ = rjo ® ^ coincides with I{t). The operator Ut defined by this integral equation is 
called again the wave operator. 

The form of the generator follows now from an operation of normal ordering. More in 
details, we start defining the flux of an observable X = X ^ 1^, where 1^ is the identity of the 
reservoir and X is an observable of the system, as jt{X) = UjXUt. Then, using the equation 
of motion for Ut and f//, we flnd that dtjt{X) = iUt[Hj'^''\t), X]Ut. In order to compute the 
mean value of this equation on the state riQ^\ so to get rid of the reservoir operators, it is 
convenient to compute flrst the commutation relations between Ut and the limiting operators 
of the reservoir. At this stage the so called time consecutive principle is used in a very heavy 
way to simplify the computation. This principle, which has been checked for many classes of 
physical models, [1], states that, if P(t) is any of these limiting operators of the reservoir, then 

[/3{t), Uf] = 0, for all t > t' . (A.5) 

Using this principle and recalling that 770 is annihilated by the limiting annihilation operators of 

the reservoir, it is now a simple exercise to compute (dtjAX)) and, by means of the equation 

''0 

{dtjt(X)) (() — {jt(L(X))) ((), to identify the form of the generator of the physical system. 

Let us now consider the case in which T > 0. In this case the state of the reservoir is no 
longer given by the vacuum ipQ. It is now convenient to use the so-called canonical representation 
of thermal states, [1]. Using the same notation of Section IV, any annihilator operator 
be written as the following linear combination 

a,,^y/Mf)c^^^ + ^4y, (A.6) 

where m{p) and n{p) are the following two-points functions, 

t 1 + Q-P^p 

m{p) = ujf^ia^jalj = ^ _ n{p) = u;p{alja^j) = j^^z^, (A.7) 
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for our bosonic reservoir, if ujp is a KMS state corresponding to an inverse temperature (5. The 
operators are assumed to satisfy the following commutation rules 

(A.8) 

while all the other commutators are trivial. Let moreover $o be the vacuum of the operators 

(a) 

Then it is immediate to check that the results in (A. 7) for the KMS state can be found, using 
these new variables, representing as the following vector state c<j^(-) = ($oi ■*^'o)- With this 
GNS-likc representation it is trivial to check that both the CCR and the two-point functions 
arc easily recovered. This representation is also called in [1] the Fock-anti Fock representation 
because of the different sign in the free time evolution of the annihilation operators ci"- and 
c^'']. Once this representation is introduced, all the same steps as for the situation with T = 
can be repeated, and the expression for the generator can be deduced using exactly the same 
strategy. 
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